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Abstract. The evolution of large-scale magnetic fields in 
disk galaxies is investigated numerically. The gasdynam- 
ical simulations in a disk perturbed by spiral or bar po- 
tential are incorporated into the kinematic calculations 
of induction equations to elucidate the effects of non- 
axisymmetric disk structure on magnetic fields. The ef- 
fects of interstellar turbulence are given as the turbulent 
diffusion of magnetic fields. The usually adopted dynamo 
mechanism of a-effect is not considered in our compu- 
tations, because it is not obvious about the actual exis- 
tence of the effect in a galaxy. Our principal concern is to 
clear how observationally and theoretically well-established 
gas flow affects the magnetic-field structure and evolution, 
without putting a lot of artificial parameters in the model. 

We have found that the density-wave streaming mo- 
tion of gas has a significant influence on the distribution 
of magnetic fields: the lines of force are well aligned with 
spiral arms due to the compressional and additional shear- 
ing flow of gas in these regions. In the inter-arm regions, 
the field lines have the finite angles with respect to the 
imposed arms, because the gaseous streamlines induced 
by spiral arms deviate from the orientation of spiral arms 
themselves. These properties of magnetic-field orientation 
across the arm are well in agreement with the results of 
radio continuum observations. All simulation models have 
resulted in the eventual decay of magnetic energy due to 
the strong turbulent diffusion. 

We have also explored the azimuthally periodic func- 
tion for the coefficient of turbulent diffusion with the max- 
imum in spiral- arm regions. This is anticipated from the 
enhancement of turbulence by young OB stars and super- 
novae. This effect of non-uniform dissipation gives rise to 
the modulation of magnetic-field structure and its time 
evolution. In particular, the dissipation rate of magnetic 
fields can be much smaller than the usually assumed rate 
of - 1/10 8 yr" 1 . 
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1. Introduction 

Measurements of linear polarization in radio continuum 
observations have revealed characteristic morphology of 
magnetic field configuration in galaxies (Beck 1993): lines 
of force are well aligned with characteristic features of a 
disk galaxy such as spiral arms and bars, and a few but 
not all edge-on galaxies prevail largely vertical fields per- 
pendicular to a disk. Attempts to derive actual directions 
of magnetic vectors are rather delicated tasks with a lot of 
ambiguities, and thus have not well succeeded to uncover 
whether field directions change along the azimuth or not in 
a galactic disk. The recent high-resolution measurements 
also allow to derive the fine structure of the large-scale 
magnetic fields in a disk: pitch angles of these fields vary 
with azimuthal angle in the systematic way (e.g. Ehle & 
Beck 1993), and the well uniform fields are most recog- 
nizable outside the optical spiral arms. Apparently, these 
observations suggest that origin of magnetic morphology 
is principally associated with a non-axisymmetric galac- 
tic structure accompanying a characteristic large-scale gas 
flow. 

It has been argued that present structure and typi- 
cal strength of galactic magnetic fields are organized by 
induction processes called galactic aw-dynamos (Krause 
et al. 1993). Galactic dynamos are invoked by the self- 
amplification of magnetic fields via helical small-scale tur- 
bulence (a-effects) and differential galactic rotation [u>- 
effects). Generally initial (seed) magnetic fields are com- 
posed of several magnetic modes, and the modes with 
finite exponential growth are supposed to predominate 
at the present epoch after some nonlinear saturation of 
growth. The conventional models are characterized by a 
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rather slow time scale of field amplification, of the or- 
der of Gyr, and attempts have been made to explain the 
present field structure, mainly based on an axisymmetric 
disk model with some additional effects of vertical flow. An 
implicit assumption in such an approach is that a galactic 
structure is only the second order effect on a magnetic- 
field structure, and dynamo experts tend to invoke the 
properties of turbulence (such as tensor a) to explain non- 
axisymmetric magnetic fields (e.g. Riidiger et al. f993). 

Basically an induction equation says that the 
magnetic-field evolution is strongly governed by the pro- 
files of velocity fields, which are typically having a time 
scale for change, of the order of f 8 yrs, in a galactic disk. 
Then it is followed that most of the present field struc- 
tures reflect the velocity fields of gas which have hap- 
pened within ~ fO 8 yrs, irrespective of slow aw-dynamo 
processes. Dynamos might give a solution for global field 
geometry (e.g. axisymmetric or bisymmetric), which is 
however not well defined from observations. It has also 
been claimed that dynamos cannot explain the strong field 
at high-redshift intergalactic clouds (Wolfe et al. 1992). 
Therefore, as most observational results suggest as men- 
tioned above, it is essential to clarify the role of a galactic 
structure in understanding the present structure of mag- 
netic fields (see Kronberg 1994 and Chiba & Lesch 1994 
for more detailed discussion about the present status). 

In the present paper, we focus on the evolution of mag- 
netic fields under the influence of non-axisymmetric dis- 
turbance, such as spirals and bars, and perform numeri- 
cal simulations to elucidate how such disturbances affect 
the evolution and structure of magnetic fields. With the 
same motivation, Panesar & Nelson (1992) have shown 
from their simulation that spiral arms strongly modify the 
field morphology. Since they have performed only one case 
due to the extensive nature of simulation, we attempt to 
consider several cases for the parameters of disturbances, 
particularly concerning the dependence on pitch angles of 
spiral arms, and spatial profiles of turbulent diffusion for 
magnetic fields. 

The plan of the paper is as follows: Sec. 2 is devoted 
to methods of simulations and model parameters. Results 
of simulation are explained in Sec. 3. Conclusions and dis- 
cussion, especially the comparison with observations, are 
given in Sec. 4. 

2. Numerical Methods and Model parameters 

2.1. Hydrodynamics in a galaxy 

In our simulations of magnetic fields in galaxies, it is cru- 
cial to obtain the realistic large-scale flow of gas, which 
is most plausible theoretically and also reminiscent of ob- 
served velocity fields of the neutral gas component as well 
as the warm ionized gas. We adopt the hydro-code of 
smoothed particle hydrodynamics (SPH), first introduced 
by Lucy (1977) and Gingold & Monaghan (1977). The 



SPH code is based on a particle method, without defining 
a computational grid. Thus it is fully Lagrangian, and is 
naturally implemented in three-dimensional calculations. 
Because of its flexible nature of the scheme, it is widely 
used for many astrophysical problems (see e.g. Monaghan 
1992; Benz 1990; Hernquist & Katz 1989). 

2.1.1. Fundamentals in SPH 

In SPH, the computation consists of a number of discrete 
fluid elements, i.e. particles, so that some local averages 
must be performed using an interpolating kernel W(r; h). 
h is the smoothing length. Then any mean value of phys- 
ical quantity, say A(r) can be defined by 



<A(r) >= J A(r')W(r-r';h)dr', 



(1) 



where the kernel W is normalized, J W(r; h)dr = 1. When 
the values of A(r) are known only at finite and discrete 
points, the integral interpolant is replaced by the following 
equation, 



<A(r) >=J2 



rn 



3=1 



A{r 3 ) 
Pj 



W(T- rj ;h), 



(2) 



where rrij is the particle mass; the density pj at each par- 
ticle's position is defined as 



N 



p = ^ j mjW(Y-r j ;h). 

3 = 1 



(3) 



In our simulation, we adopt the kernel based on spline 
functions (Monaghan & Lattanzio 1985), 

x (l - (3/2)a; 2 + (3/4)a; 3 for < x < 1 
W(r;h) = — { (l/4)(2-a;) 3 for 1 < x < 2 (4) 



TTh 6 







for x > 2, 



where x — r/h. Here the smoothing length is variable in 
each particle, hj, according to the density at the particle's 
position (see below). 

The hydrodynamical equations in our model are given 

by 



dYi 



(5) 



dt 

^p- = -— VP 4 + q 4 - V$ - 2fl p x v - n p x fl p x r,(6) 

dhi hi dpi 

dt 3pi dt ' 



(7) 



where Pj is the pressure, $ the gravitational potential, 
and q, is the artificial viscosity to treat the formation 
of shock waves. Since our aim is to obtain the charac- 
teristic structure of velocity fields in a galaxy, where the 
effect of self-gravity of gas is minor (except at the galac- 
tic center), we neglect the contribution of gas itself in $ 
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but consider the external galactic potential provided by 
galaxy mass as a whole. Q p denotes the pattern speed of 
non-axisymmetric disturbance in a galactic structure (see 
next section). Therefore we simulate the gas dynamics in 
a frame rotating with f2 p , where such disturbance is ob- 
served to be stationary. 

The terms of pressure gradient and artificial viscosity 
are written as: 



Pi 



N P 
"i 



3 = 1 



i 



x- [ViWin^hi) + ViWir&hj)] , (8) 

where Uy denotes the contribution of artificial viscosity 
given by 



n, 



and 



My 



■ farv -r tf <0 







Pij 



for Vy ■ Tij > 0, 



*l + v 2 



(9) 



(10) 



Here, — (hi + hj)/2, p t j = (p l + pj)/2, and the 
average sound velocity, = (cj + cj)/2. In the ex- 
pression for Ily , eq.(^), the first term corresponds to a 
bulk viscosity while the second term is equivalent to the 
Von Neumann-Richtmyer viscosity. For the parameters 
(a,(3,r]), we adopt, a = 1, /3 = 2, and rj 2 — O.Ol/i^-, 
which reproduce reasonably some test calculations such as 
one-dimensional shock-tube problem and two-dimensional 
wind-tunnel problem. The factor / is introduced to pre- 
vent the excessive shear viscosity that appears in the stan- 
dard form of viscosity given in equs.(^) and ([To|). In our 
computation, we use the form for / described by Benz 
(1990) as, 



/ = 



< V • v > 



< V • v > 



< V x v > 



(11) 



•O.OOlcy/Zly 

which vanishes if the quantity < V x v >, the curl of 
velocity is large, suppressing the shear viscosity. 

The interstellar gas is modeled in terms of the isother- 
mal equations of state, P = c 2 p. 

The basic equations eq.(5), eq.(6), and eq.(7) are inte- 
grated in time, using the leapfrog method. The time step 
At is controlled by the Courant condition: 

hi 

I , (12) 

hi\V ■ v| + Cj + 1.2(acj + f3maxj\fiij\) 



At = 0.3mm 



2.1.2. Galaxy models 

The axisymmetric part of the gravitational potential in 
a disk plane is represented as the Toomre disk (Toomre 
1963) : 



0o(r) = (r 2 + a 2 ) 

a 



-1/2 



(13) 



which gives the maximum rotational velocity v ma x — 
(4/27) 1 / 4 c/a at a radius r max = ay/2. For the vertical 
dependence of the potential, we adopt a parabolic form, 
vz 2 /2, where v is the vertical frequency for motions of 
stars. We use the characteristic value in the solar neigh- 
borhood, v — 3.2 x 10" 15 s" 1 (Binney & Tremaine 1987), 
for all radii. Of course, this value as well as the parabolic 
approximation does not hold all over a disk, especially at 
high z where this approximation is invalid. The gas flow 
derived from this galactic potential is however well repre- 
sentative of actual flow for the present purpose to simulate 
the magnetic-field evolution. 

Initially the gas disk has the approximately uniform 
density with a finite thickness as determined by the verti- 
cal dependence of the potential given above and the sound 
speed of gas, c, which is assumed to be constant, 8kms _1 . 
The disk is set up in centrifugal equilibrium. The param- 
eter a in eq. @ is 2.6kpc in our simulation, and c is de- 
termined so that v max = 200 kms _1 kpc _1 . 

Non-axisymmetric disturbances are introduced to a 
disk in the first 1.6 x 10 s yrs of time evolution, with the 
form, 



<f>l(r, <p) = (j)o(r)e(r) cos [2 ln(r/r ) cot ?A + 2<p], 



(14) 



and are fixed hereafter (Matsuda & Isaka 1980; Johns & 
Nelson 1986). ip is the pitch angle of spiral disturbance 
with a logarithmic two-armed form, and ro, arbitrary pa- 
rameter, determines the phase of spiral. e(r) is represented 



e(r) = eo-s(l + s) 



-3/2 



(15) 



which reaches the peak value Sq^JaJTJ at a radius r max . 
In the simulation for a bar perturbation, the term propor- 
tional to cot-0 in eq.(p^) is omitted. 

Figure 1 shows the rotation curve and some frequency 
distribution in our model. It follows that the strength of 
velocity shear is strongest at r ~ 3kpc. 

In order to obtain the characteristic velocity fields for 
solving the magnetic-field evolution, the computations of 
gas dynamics are halted after a few hundreds Myr when 
the gas flow becomes approximately steady, and then the 
velocity fields are evaluated as will be explained in the 
next section. We note that a small-scale time-dependent 
modulation in gas dynamics is observed (cf. Johns & Nel- 
son 1986), but the obtained velocity profiles by the above 
procedure are enough characteristic to see their effects on 
magnetic fields. As an example, the velocity-field vectors 
for the pitch angle 20° is presented in Fig. 2. The profiles 
of gas density for the spiral pitch angles, 10° and 20°, and 
the bar perturbation are shown in Fig. 3. The solid lines 
mark the minima of the gravitational potential of spirals. 
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Fig. 1. The rotation curve (a) and frequency distribution (b) 
in the galaxy model 

2.2. Solving induction equations 

The kinematic evolution of the magnetic field is analyzed 
using the time dependent solution of the field transport 
equation with diffusion (the induction equation), 

dB/dt = V x (v x B) - V x rj (V x B), (16) 

where r\ is the diffusion coefficient, B and v are the mag- 
netic and velocity fields, respectively. The assumed ap- 
proach allows us to investigate qualitatively the effects of 
spatially variable r\ on the evolution of magnetic fields in 
a disk. However, the more precise solution of this problem 
is also possible (Kichatinov & Riidiger 1992). 

The problem is solved numerically in 
three-dimensional (3D) rectangular grid of points, where 
the XY plane is the galactic plane and the Z-axis is the 
axis of galactic rotation pointing to the North. The rect- 
angular size in X and Y directions is 20 kpc, where the 
grid interval is 400 pc. Along the Z axis the box has 3 or 
4 kpc depending on the assumptions concerning the dif- 
fusion coefficient (see below) and with the spatial step of 
250 pc interval. The time interval is At = 0.22 • 10 6 yrs. 
The fourth order explicit Runge-Kutta method is used for 
the discretization in time. The space derivatives are ap- 
proximated by the spline functions of the third order. As 
the boundary conditions, the continuity of the third order 
derivatives of the spline functions has been used. Physi- 
cally it means that impenetrable walls have been applied. 



Fig. 2. The velocity field vectors in the case of the spiral pitch 
angle V = 20° 

This fact causes the necessity of the velocity field decreas- 
ing near the edges. 

The prevalent number of numerical methods applied 
to the induction equation do not satisfy the divergence- 
free condition of the magnetic field. There are several ap- 
proaches which can be used to solve this difficult problem 
(Elstner et al. 1990; Evans & Hawley 1988; Schmidt- Voigt 
1989). We have adopted the "divergence-cleaning" method 
elaborated by Schmidt- Voigt (1989) which is suitable in 
our case (see discussion in Schmidt- Voigt 1989 and Evans 
& Hawley 1988). Following Schmidt- Voigt, in order to di- 
minish the errors of the non-zero magnetic divergency, we 
introduce a potential <J>f, which is a solution of the Poisson 
equation, 

V 2 * b (x,y,z,t) = -(V • B). (17) 

After each time step, this equation is solved using the 
Gauss-Seidel iteration, and we replace the magnetic field 
Bby 

B' = B + V $ 6 . (18) 
The iteration ends when the following condition is fulfilled: 

|(V-B)| • | Ax | / |B| < 0.01. (19) 

This is what is necessary to avoid the unphysical solutions 
(Schmidt- Voigt 1987). 

To incorporate the particle-based flow in SPH into 
the grid-based code for magnetic-field evolution described 
above, the averaging of velocity fields in grids is performed 
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Fig. 3. The profiles of gas density (grayplot) in the cases of the 
spiral pitch angle ip — 10° (a), if) = 20° (b), and the bar (c). 
The solid lines denote the minima of the gravitational potential 
of spirals 



using the smoothing length hi in SPH and Gaussian kernel 
W g (Habe, et al. 1991): 

1 N 

Vyfc = m n v n Wg(Anj k ,n; h n ), (20) 

Pijk , 

where Ar ijkin = \(x n - x;) 2 + (y n - yj) 2 + (z n - 
Zfe) 2 ] 1 / 2 , pij k = J2n m n W g( Ar vk,n;h n ), and W g = 
l/ir 3 / 2 /hn exp[— (r/h n ) 2 ]. Here the indexes i,j, k are for 
grids, and n for particles. The reason to adopt the Gaus- 
sian kernel in this procedure is that due to the finite num- 
ber of particles (N — 2000), there would be the lack of ve- 
locities at some grids in the spline case, because W defined 
in eq.(|]) is zero for r/h > 2. However, the different choice 
of kernel for velocity fields at grids does not introduce the 
fundamental change in the characteristic properties of gas 
flow obtained in SPH. 

The computations have been performed using Cray Y- 
MP in HLRZ Jiilich and SPARC 10 in the Max-Planck- 
Institut fur Radioastronomie in Bonn. 

2.3. Model parameters 

Having combined 3-D SPH code with 3-D induction code, 
we simulate the evolution of magnetic-field structure and 
strength in a model spiral galaxy. Although there are 
many factors to affect velocity- field properties, owing to 
the limitation of available computation time, we will focus 
only upon some characteristic gas flows. 

The model parameters adopted in the computations 
are summarized in Table 1. Our models are divided into 
2 groups depending on the initial condition for the mag- 
netic field: axisymmetric and uniform field configurations. 
Though more general would be the superposition of sev- 
eral magnetic modes (Panesar & Nelson 1992), these ini- 
tial conditions are enough to see the possible effects of 
non- axisymmetric velocity disturbances, and also the pre- 
dominant field configurations revealed in radio observa- 
tions. In Table 1, the models Ai(i = 1 — 4) are devoted 
to the initially axisymmetric magnetic-field configuration. 
Models Bi(i — 1 — 7) correspond to initially uniform (par- 
allel to the Y axis) magnetic-field structure. The inten- 
sity of the field is 1 pG in both cases from the reason 
that the microgauss order of initial magnetic fields after a 
disk formation has been suggested by the observations of 
damped Lya clouds (e.g. Kronberg et al. 1992, Wolfe et al. 
1992), and an idea to explain such observations was given 
in Lesch and Chiba (1994). Note that since the present 
concern is devoted to linear processes and the resultant 
structure of magnetic fields, the actual value of initial 
field strength does not change our results. The vertical 
velocity of gas is forced to be zero in all experiments for 
the sake of simplicity. For the diffusion coefficient 77, we 
adopt the magnitude resulting from the small-scale tur- 
bulent gas motions. Its basic value in a galactic plane is 
770 = 5.7 • 10 27 cm 2 s . This rather high value for 770 
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Table 1. The model parameters 
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compared to the usually adopted value is to reduce the 
computational time for solving the induction equation. We 
have confirmed from the simulation with smaller 770 that 
the essential results do not change at all. We note that 
this order of magnitude for 770 is close to the upper limit 
of the turbulent diffusion coefficient for the interstellar gas 
suggested in the literature (Lesch 1993). 

We have introduced the z-dependent diffusion coeffi- 
cient, rj(z): the presence of the active corona in the galaxy 
with high temperature (10 6 K) and the high turbulent 
velocity causes the increase of rj in the direction perpen- 
dicular to the galactic plane (Camenzind & Lesch 1994; 
Spencer & Cram 1992). From this reason we have applied 
the following form for 77 (Spencer & Cram 1992): 

{770 for z < H dif 

T] ■ exp ((z - H dif )/h) for H dit < z < iJ max (21) 
for Z > ffmax, 

where 770 is our basic value of the diffusion coefficient, 
-^dif ( = 0-4 kpc) is the scale height of the interstellar gas, 
and h is the scale length of the increase of 77. In our case, 
h = 500 pc, so that 77 increases about 2.5 times that 
in the galactic plane. In order to stabilize the code near 
the boundaries in the vertical direction, the diffusion co- 
efficient above H max = 1 kpc is constant with a value 
?7max — 1.4- 10 28 cm 2 s _1 . The average value of 77 adopted 
above suggests a magnetic Reynolds number R m = 55 
[defined as R. m — v max L/rj, where i> ma x is the maximum 
value of the velocity field (200 km s _1 ), and L the radius 
of the model galaxy (10 kpc)]. Relatively high value of 77 
causes the characteristic evolution time of magnetic field 
less than 10 9 yrs, thereby allowing us to reduce the ac- 
tual time of computation. However, for one experiment, 
we have adopted the diffusion coefficient ten times smaller 
than the basic one, for which we have applied At = 10 5 
yrs and assumed the uniform magnetic field as an initial 



field configuration. As a result, we obtain R m = 430 (Ta- 
ble 1, B7). 

In order to take into account the influence of spiral- 
arm activity on the diffusion coefficient, we have adopted 
the following form for rj(r, (p): 

V(r,tp) - 7/0 ~ 2 T?mm [1 + cos{2y + 2 ln(r/r ) cot 7^}] 
+77 

mm ? 

(22) 

where r\ m im the minimum value of rj, is equal to 770 /10 
in our model. Equation (22) is equivalent to eq.(14), i.e. 
we assume the two-armed logarithmic spirals for 77. Since 
the turbulent velocity of interstellar gas, < v 2 > 1 / 2 , would 
be high in arm regions due to efficient energy input from 
young OB stars and supernovae, the turbulent diffusion, 
77 ^< v 2 > 1 / 2 / where I is the correlation length of turbu- 
lence, would be also large (Ko & Parker 1989). Therefore, 
we assume in our simulations that 77 is maximum in arm 
regions. Of course, it is not obvious as to whether this as- 
sumption is realistic or not. However, the depolarization 
of radio continuum emission preferentially in arm regions 
may support this idea (e.g. Horellou et al. 1993; Neiningcr 
et al. 1993). The magnetic Reynolds number is about 178 
in this case (B5). 

Two experiments (A2 and B6) have been made for the 
diffusion coefficient variable in the galactic plane as well 
as in in the z direction accordingly to the above formulaes 
(21 and 22). The averaged value of 77 gives as a result 
R m = 145. 

To avoid numerical instabilities related to the steep 
diffusion gradients, rj has a constant value in the central 
part of the disk (to the radius of 2 kpc) in all cases. 

3. Results 
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Fig. 4. Intensitites of magnetic fields (grayplot) and magnetic vectors for the model Al at t = 1.1 • 10 8 yrs (a) and at t = 2.2 • 10 s 
yrs (b), and for the model A2 at the same epochs (c) and (d). The solid lines denote the minima of the gravitational potential 
of spiral arms 



3.1. The 3D evolution of magnetic field: the case of ini- 
tially an axisymmetric magnetic field 

The calculations with initially a ring magnetic field struc- 
ture have been done for three types of gaseous velocity 
fields (Table 1; A): two- armed spiral disturbances with the 
pitch angle ip being 10° and 20°, and the bar perturba- 
tion. By these experiments, we attempt to clarify the net 
influence of non-axisymmctric disturbance on an axisym- 
metric magnetic-field distribution. In all cases, the mag- 
netic field follows closely the gas motions, thus prevailing 
similar structures to spirals or a bar at early evolution- 
ary stages. As time progresses, the field is expelled out 
from the inner part of the disk due to the adopted strong 
turbulent diffusion. 

Fig. 4 presents the grayplots of the intensity and the 
vectors of the magnetic field in the galactic plane, for the 
models Al with rj{z) (a and 6), and A2 with r](r,(f>,z) 
(c and d). Fig. 4a and 4c corresponds to the evolution- 



ary stage of t = 1.1 • 10 8 yrs, while Fig. 4b and 4d for 
t = 2.2 • 10 8 yrs. The solid lines denote the minima of the 
gravitational potential of spirals. For both experiments, 
the initially circular field is transformed into the spiral 
structure (a and c) and then diminishs by the strong dif- 
fusion (6 and d). The intensity maxima of the magnetic 
field are visible slightly outside the potential minima (i.e. 
down-stream side), which are close to the regions where 
the velocity shear is the highest one. Comparing Fig. 4 
with the gas density map for ip = 10° (Fig. 3a), it fol- 
lows that the maxima of the magnetic field coincide with 
those of gas density, though in the magnetic field case, the 
outer magnetic layers of the galactic disk are much more 
dispersed. 

Fig. 4 also shows the visible difference between the ex- 
periment with the diffusion coefficient r\ (z) and r\ (r, <f>, z): 
the applied contrast of r\ between the interarm and arm re- 
gions causes higher magnetic field amplification near and 
inside the spiral arms. It means that the variability of the 
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Fig. 5. The grayplots of pitch angle distributions of magnetic 
fields for ip — 10°. The dark (white) areas denote the positive 
(negative) pitch angles (the colors range from -25° to 25° every 
2°). The solid lines mark the gravitational potential minimum, 
(a): Al with r)(z). (b): A2 with ij(r,<p,z). Fig. 5(c) shows the 
pitch-angles of velocity fields 



diffusion coefficient can play an important role in the evo- 
lution of galactic magnetic fields. 

One of the most important goals of our simulations 
is to answer the question how closely the magnetic field 
follows the interstellar gas motions. In order to analyze 
it, we introduce the pitch angles of magnetic and veloc- 
ity fields, defined as the angle between the field vector and 
the azimuthal direction in the galactic plane. We note that 
the obtained distribution of magnetic pitch angles is not 
directly compared with observational results, because the 
model does not take into account some observational ef- 
fects such as the beam smoothing or the Faraday rotation. 
Our aim is to clear the physical connection between mag- 
netic and velocity fields, which is free from ambiguities in 
interpretation of observations, and tremendous number of 
freedom in parameters depending on galaxy inclinations, 
disk thickness, effects of small-scale magnetic fields, and 
so on. 

Fig. 5 shows the grayplots of the magnetic pitch an- 
gles at evolutionary time t = 1.1 ■ 10 8 yrs in the cases of 
r] (z) (Al) (Fig. 5a) and r\ (r,(p,z) (A2) (Fig. 5b), as well 
as, the pitch-angle distribution of the velocity field with 
ip = 10° (Fig. 5c). The dark areas correspond to posi- 
tive pitch angles - the vectors are directed outward from 
the disk compared to the azimuthal direction, while the 
white ones for negative pitch angles - the vectors are di- 
rected inward. The absolute maximum values of the mod- 
elled pitch angles are below 30° what is in good agreement 
with observational data (e.g. Neininger 1992, Beck 1993). 
The solid line, as in Fig. 4, shows the minimum of gravi- 
tational potential of spirals. It is easily seen that all three 
maps for pitch angles are very similar: the velocity field 
as well as the magnetic one is directed outward in front 
of the potential minimum (i.e. upper-stream side), then is 
aligned along the spiral arms, and finally is directed in- 
ward the galaxy. The highest contrast in the intensities 
of pitch angle is obtained for i] (r,<fi,z), manifesting the 
greatest evolutionary changes of magnetic-field structure 
due to the gas motions and the diffusion. 

The numerical experiment using spiral velocity fields 
with ip = 20° (model A3) has resulted in the higher 
amplification of the magnetic-field intensity (Fig. 6a and 
6b, the legend is the same as in Fig. 4a and 4b). The 
maxima are distributed just outside the potential min- 
ima, more significantly than in the previous case. As ob- 
tained for ip — 10°, the magnetic maxima coincide with 
the gas density maxima (see Fig. 3b). Fig. 6a is drawn at 
t = 1.1 • 10 8 yrs. When we apply two times longer evo- 
lution (Fig. 6b), the magnetic flux decreases more slowly 
than in the Al or A2 computations. Fig. 6c and 6d present 
the pitch angle grayplots for the magnetic and velocity 
fields, respectively. The description is the same as for 
Fig. 5. The distribution of both fields is almost identical. 
The magnetic as well as the velocity field are directed out- 
ward from the disk in front of the spiral potential, then are 
tightly aligned with the arms, and finally they point in- 
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Fig. 6. Results of the model A3, (a) and (b) show the magnetic-field intensity (grayplot) and vectors at t — 1.1 • 10 8 yrs and 
t — 2.2 • 10 8 yrs, respectively, (c) and (d) show the pitch-angle distributions of magnetic and velocity fields, respectively. The 
solid lines denote the gravitational potential minimum 



ward. Thus the situation is similar to the i\) = 10° case, 
but now, due to the higher non-axisymmctricity in the 
gas flow, the shear is stronger, causing higher magnetic 
field amplification. However the absolute maximum value 
is again less than 30°. 

The last case for initially a ring magnetic field is the ex- 
periment with the bar perturbation on the gaseous veloc- 
ity field (model A4; Fig. 7). The magnetic intensity max- 
ima are not large (Fig. 7a) and they are fully dispersed 
at the evolutionary stage of t = 2.2 • 10 8 yrs (Fig. 7b). 
In the central part of the disk, the strong diffusion sweeps 
the flux quite significantly. The shear in the gas velocity is 



not high so that the amplification of magnetic field is low. 
The distribution of the magnetic pitch angles, which are 
also similar to those of the velocity fields, proves that the 
magnetic field is aligned along the bar and arms, changing 
its direction from outward to inward across the bar and 
spiral structures. 

The magnetic energy density decreases in all cases pre- 
sented, reaching the lowest value after 5 • 10 s yrs due to 
the strong turbulent diffusion. 
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3. 2. Results for initially a uniform magnetic fields 

The calculations of magnetic field evolution with initially 
a uniform (parallel to the Y axis) magnetic field have been 
performed for four velocity field configurations (Table 1; 
models Bi): the pure circular motion with differential ro- 
tation, spiral disturbances with ip being 10° and 20°, and 
the bar perturbation. 

For all experiments, the magnetic field is wound into 
the bisymmetric structure, causing the initial growth of 
the magnetic energy density (Fig. 8). As time progresses, 
the flux in the inner part of the disk is quickly swept 
out by the turbulent diffusion, and the magnetic energy 
decreases. Fig. 8 demonstrates the time evolution of the 
magnetic energy density e normalized to the initial en- 
ergy density eo- The curves are presented for four chosen 
cases with Reynolds number R m = 55 (Bi(i = 1 — 4)), 
R m = 178 (55), R m = 145 (B6) and the case with 
R m = 430(57). On the contrary to the experiments with 
initially a ring magnetic field, the smallest growth of e oc- 
curs for the spiral velocity with ip = 10°, while the largest 
one for the bar disturbance (the cases Bi(i = 1 — 4)). It 
follows that for the flow with higher non-axisymmetricity, 
the winding and dispersing of the magnetic field is faster. 
The case with r](r, <fr, z) and ip = 10° (BQ) gives the high- 
est increase of the magnetic energy, mainly due to the 
smaller average value of the diffusion coefficient. The ex- 
periment B5 with the diffusion coefficient variable in the 
galactic plane, but constant in the z direction (what could 
have place for galaxies without the active corona) reveals 
two times bigger growth of e and farther maximum with 
slowly decreasing slope manifesting sustaining of the mag- 
netic field longer than the computation time. 

The calculations for the largest Reynolds number of 
R m = 430 (B7) are the most realistic to the real galaxy 
conditions. The growth of the magnetic energy density is 
about three times greater than for the same experiment 
with R m — 55 (B2). The curve in this case does not de- 
crease smoothly but shows two farther maxima. This fact 
is well demonstrated in the maps of the magnetic-field 
vectors shown in Fig. 9. After 10 8 yrs, the clear bisym- 
metric structure is visible with the slightly dispersed in 
the center (Fig. 9a). This configuration, though more dif- 
fusive, can be still observed at t = 2 • 10 8 yrs (Fig. 9b) 
and t — 3.5 • 10 8 yrs (Fig. 9c). In Fig. 9d, the magnetic 
flux is completely swept from the central part of the disk, 
whereas there is still a kind of ring-like field structure in 
the regions near the edges. It means that for experiments 
with smaller diffusion coefficient, magnetic field structures 
will be longer affected by the non-axisymmetric gas flows. 



Fig. 7. Distributions of magnetic intensity and vectors in a bar 
perturbation (model AA) at t = 1.1-10 8 yrs (a) and t = 2.2-10 8 
yrs (b) 



4. Conclusions and discussion 

The evolution of large-scale magnetic field in the galac- 
tic disk affected by non-axisymmetric structures like spi- 
rals and bars has been demonstrated using the 3D nu- 



Fig. 8. Time evolution of magnetic energy density e normalized 
to the initial energy density eo for the models Bi. Circ, bar, p20 
and plO denote the velocity field structures: circular, with bar 
perturbation, with the pitch angle of 20° and 10°, respectively 
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Fig. 9. Distributions of magnetic vectors for the model B7 at 4 evolutionary epochs 



mcrical simulations. First, to obtain the realistic velocity 
fields for calculations of magnetic fields, the hydrodynam- 
ical simulations using SPH scheme have been performed 
for three potential disturbances in an axisymmetric disk: 
tightly wound spiral with pitch angle ip of 10°, open spiral 
with ip = 20° and bar perturbation. Then, the resultant 
velocity fields of gas at ~ 1.6 ■ 10 8 yrs have been used 
to simulate the evolution of large-scale galactic magnetic 
fields. Our calculations of magnetic fields are devoted to 
a disk structure which is well established observationally 
and theoretically: some other effects such as vertical gas 
flows are neglected because these are tediously increas- 
ing parameters for models. Two basic configurations for 
the initial magnetic field have been chosen: axisymmetric 



(Table 1, Ai) and uniform (parallel to the Y axis) (Ta- 
ble 1, Bi). The adopted field configurations are enough 
to demonstrate the net effects of non-axisymmetric flow 
perturbations on magnetic field structure. 

4-1. Magnetic field structure 

The robust conclusion from our simulations is that the 
magnetic field follows closely the interstellar gas motions; 

The comparison between magnetic energy and gas den- 
sity in the case of initially ring magnetic field (Ai) shows 
that the maxima of magnetic energy coincide with those 
of gas density (Figs. 3, 4, and 6 a and b), and this is also 
the case for initially uniform field. Both maxima are lo- 
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eating slightly outside the minima of spiral potential (i.e. 
down-stream side) within the corotation radius of 9.1 kpc, 
being the characteristic non-linear response of interstellar 
gas (if we consider the effects of magnetic pressure con- 
sistently, the region of density maxima can be broadened, 
see e.g. Tubbs 1980). 

The pitch angle distributions of both velocity and mag- 
netic fields also indicate the tight correlation between 
the interstellar gas motion and magnetic-field structure 
(Fig. 5a, b and c; Fig 6c and d): both distributions of pitch 
angles are very similar even for the longest evolutionary 
time. This suggests that the characteristic gas motions in- 
volved in a non-axisymmetric disk are dominant enough to 
organize the field structure even in the presence of turbu- 
lent dissipation. The actual spatial variation of field lines 
of force is summarized as follows: In front of the spiral 
arm (i.e. upper-stream side), the orientation of magnetic 
field lines is slightly deviating from the arm, with the ab- 
solute value of pitch angle ip being small or positive ip 
(i.e. more leading sense than the trailing spiral arm it- 
self). In the arm, the field lines are well aligned with the 
arm, i.e. the absolute value of ip is maximum with the neg- 
ative sign. Then, outside the arm, the orientation of field 
lines is changing gradually from the trailing sense in the 
arm to the leading sense as the gas stream leaves the arm, 
so that the field crosses with the next arm again with a 
finite angle. Thus, the variation of the magnetic pitch an- 
gles with azimuth is doubly-periodical, and the properties 
of the pitch-angle variation are essentially in agreement 
with the observational results for NGC6964 (Ehle & Beck 
1993), M51 (Neininger 1992), M83 (Neininger et al. 1991), 
and M81 (M.Krause, private communication). We remark 
that the definition of pitch angles by radio observers is dif- 
ferent from the usual one, and thus one should be careful 
when saying the sign of pitch angles. Also, our estimate of 
the field orientation does not include some observational 
effects such as the beam smoothing, Faraday rotation and 
the effects of small-scale turbulent fields, so that the ob- 
served distribution of the field lines obtained above may 
be somewhat different but the general properties of field 
structure across the arm may be unchanged. 

We have clarified that the actual structure of mag- 
netic field is well correlated with the optical structure of a 
galactic disk as radio observations have revealed. We are 
thus puzzled how other approaches based on conventional 
acj-dynamos with only an axisymmetric disk can explain 
such universal correlation without relying on spiral-arm 
disturbance. 

4-2. Time evolution of magnetic field energy 

All of our simulated models result in the decay of mag- 
netic energy density in the absence of field generation 
mechanisms. In the cases of initially ring fields (Ai), the 
magnetic energy density, e, simply decreases with the evo- 
lutionary time due to the weakness of the coupling be- 



tween convectional motion of gas and magnetic fields, 
|v x B| <C I • As mentioned before, each model suggests 
the amplification of magnetic field in the arm due to com- 
pressional and additional shearing motions of gas induced 
by the arm; the strongest amplification is obtained for 
the spiral potential when the pitch angle tp of 20° (^43). 
On the other hand, the experiments Bi have resulted in 
the wound magnetic-field configurations causing the ini- 
tial growth of the magnetic energy density e in the first 
10 8 yrs of evolution (Fig. 8). The fast increase of e is fol- 
lowed by quick decaying of its value due to the winding of 
magnetic lines and the resultant stimulated effects of dif- 
fusion. The higher non-axisymmetric gas flow causes the 
lower growth of the magnetic energy density because of 
the stronger reconnection of the magnetic field lines. In 
two given cases (B5 and B7), the longer sustaining of e 
is observed: in the case B5 with the diffusion coefficient 
variable in the galactic plane and in the experiment B7 
with the high value of the magnetic Reynolds number of 
430. The case B5 has the smaller average value of r\ and 
higher R m (178) than others, so the slower decrease of e 
is understandable. 

The anti-dynamo theorem says that any two- 
dimensional flows result in the decay of magnetic fields 
(e.g. Zeldovich et al. 1983). This is due to the absence of 
the coupling between the field components along the plane 
of flows and the perpendicular field component. No cou- 
pling is eventually followed by the catastrophic reduction 
of the characteristic scale of magnetic fields, and thus the 
dissipation. As is well known, the so-called dynamo cycle 
relying on helical motion of turbulence can introduce the 
coupling among three components of magnetic fields. A 
disadvantage in galactic dynamo theory is that it is not ob- 
vious concerning the actual existence of helical interstellar 
turbulence from both theory and observation (cf. Ferriere 
1992; Hanasz & Lesch 1993; Kaisig et al. 1993). Spencer & 
Cram (1992) have made the important progress by propos- 
ing a new field-generation mechanism relying on global 
galactic winds. This approach is however also subject to 
the uncertainty as to whether such preferential winds can 
be realized in an actual galaxy. It is also worthwhile to 
note that their model is based on an axisymmetric disk, 
and thus suggesting the generation of non-axisymmetric 
magnetic fields in order not to break the anti-dynamo the- 
orem. The emergence of non-axisymmetric fields will give 
rise to the existence of non-axisymmetric coupling terms 
between velocity and magnetic fields, which will eventu- 
ally affect the basic mechanism of field generation from 
an axisymmetric disk. The general analysis of such a sys- 
tem is now underway (Spencer, private communication). 
On the other hand, observations of edge-on galaxies sug- 
gest the mass outflow from a disk to a halo (e.g. Bloemen 
1991), though it is not clear as to whether such vertical 
flow is occurring steadily or intermittently. 

Therefore, it is necessary to consider the general prob- 
lem of magnetic-field evolution containing the realistic ver- 



K. Otmianowska-Mazur & M. Chiba: Galactic magnetic fields 



13 



tical gas flow and the non-axisymmetric velocity distur- 
bance. If the condition for a dynamo is fulfilled in such 
a system, the growth of magnetic fields can be very fast 
when the efficient amplification due to non-axisymmetric 
disturbance, e.g. the wave coupling between density and 
magnetic waves, is considered (Chiba & Lesch 1994), 
which we have neglected here by suppressing the propaga- 
tion effect of spiral arms with azimuth in the simulations. 

4-3. Effects of azimuthally varied rj 

In the present numerical experiments, we have explored 
the possibility that the dissipation effects for magnetic 
fields are spatially varied in a galactic plane. Usually the 
coefficient rj of turbulent diffusion has been assumed to be 
uniform in a plane for galactic dynamos. Ko & Parker 
(1989) have examined the time-dependent dynamo on 
the assumption of spatial variability in dynamo efficiency 
driven by young OB stars and supernovae in spiral arms. 
The energy transformation from stellar winds and explo- 
sions to interstellar medium gives rise to high velocity dis- 
persion of turbulence, thereby enhancing the dissipation 
of magnetic fields. Mestel & Subramanian (1991) consid- 
ered the enhancement of helical turbulence (a-effects) in 
spiral arms. Though the exact expression for r)(ip) is un- 
known, it is interesting to see how the effect is non-trivial 
for magnetic-field evolution. 

The model with the diffusion coefficient variable in the 
galactic plane (A2) has resulted in higher amplification of 
magnetic intensity in the spiral arms than in the case for 
rj constant in the plane (Al). This is well explained by the 
transport of magnetic flux from the region where diffusion 
is high (in our case, the minima of potential perturbation) 
to the region where diffusion is weak. The gas streaming 
motion across the arm prohibits the leakage of magnetic 
flux towards the upper-stream side, so the magnetic inten- 
sity is magnified at the down-stream side of spiral arms. If 
we consider the generation of small-scale magnetic fields 
from such intense turbulence, the total strength of mag- 
netic fields can be the highest in the arm regions. Further- 
more, the contrast in magnetic pitch angles is also magni- 
fied in this case (Fig. 5b), probably due to the higher cou- 
pling between velocity and magnetic fields in the inter-arm 
regions where rj is assumed to be small. 

One more noticeable point for the case of spatially var- 
ied rj is that the decay of magnetic field energy is made 
much slower (Fig. 8). Therefore, more realistic may be the 
rather slow dissipation of large-scale magnetic fields com- 
pared to the usually voiced value of dissipation time scale 
of ~ 10 s yr. Further studies about this point are necessary 
before drawing such conclusions. 
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